%--------------------------------------------------------------------------
%
% Ephemeris computation using variable-order Radau IIA integrator with
% step-size control
%
% Last modified:   2018/02/11   M. Mahooti
%
%--------------------------------------------------------------------------
function Eph = Ephemeris_new(Y0, t0,tf)

options = rdpset('RelTol',1e-13,'AbsTol',1e-16);
[t,yout] = radau(@Accel,(t0:60:tf),Y0,options);
Eph(:,1) = t;
Eph(:,2:7) = yout;

